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Long-duration gamma-ray bursts (GRBs) are an extremely rare outcome of the collapse of 
massive stars and are typically found in the distant universe. Because of its intrinsic luminosity 
(Z. ~ 3 x 10 53 ergs per second) and its relative proximity (z = 0.34), GRB 130427A reached the 
highest fluence observed in the y-ray band. Here, we present a comprehensive multiwavelength 
view of GRB 130427A with Swift, the 2-meter Liverpool and Faulkes telescopes, and by other 
ground-based facilities, highlighting the evolution of the burst emission from the prompt to the 
afterglow phase. The properties of GRB 130427A are similar to those of the most luminous, 
high-redshift GRBs, suggesting that a common central engine is responsible for producing 
GRBs in both the contemporary and the early universe and over the full range of GRB 
isotropic energies. 


G amma-ray burst (GRB) 130427A was the 
brightest burst detected with Swift (7), as 
well as with several y-ray detectors on- 
board other space missions. It was also the brightest 
and longest burst detected above 100 MeV, with 
the most energetic photon detected at 95 GeV (2). 
It was detected with Fermi-Gamma-ray Burst Mon- 
itor (GBM) (3) at T 0 GBM = 47:06.42 UT on 27 April 
2013. Hereafter, this time will be our reference time 
T 0 . The Burst Alert Telescope (BAT) (4) onboard 
Swift triggered on GRB 13 0427 A at time (t) = 
51.1 s, when Swift completed a preplanned slew. 
The Swift slew to the source started at t = 148 s 
and ended at t = 192 s. The Swift Ultraviolet Op- 
tical Telescope (UVOT) (5) began observations at 
t = 181 s, whereas observations with the Swift 
X-ray Telescope (XRT) (6) started at t = 195 s (7). 
The structure of the y-ray light curve revealed by 
the Swift-BAT in the 15- to 350-keV band (Fig. 1) 
can be divided in three main episodes: an initial 
peak, beginning at t = 0. 1 s and peaking at t = 
0.5 s; a second large peak, showing a complex 
structure with a duration of ~20 s; and a third, 
much weaker episode starting at t ~ 120 s, show- 
ing a fast rise/exponential decay behavior. The 
overall duration of the prompt emission was ^90 
(15 to 150 keV) = 276 ± 5 s (the time containing 
90% of the fluence) calculated over the first 1830 s 
of BAT observation from T 0 GBM . During the 
early phases of the y-ray emission, strong spectral 
variability was observed (Fig. 1). A marked spec- 
tral hardening was observed during the prompt 
main event. With a total fluence F = (4.985 ± 
0.002) x 10 4 erg cm -2 in the 15- to 150-keV band, 
GRB 13 0427 A reached the highest fluence ob- 


served for a GRB by Swift. The 0.02- to 10-MeV 
fluence measured by Konus-Wind (5) for the main 
emission episode (0 to 18.7 s) is (2.68 +/- 0.01) x 
1CF 3 erg cm -2 , with a spectrum peaking at E ^ k 
1 028 ± 8 keV, whereas the fluence of the emission 
episode at (120 to 250 s) is ~9 xio -5 erg cm -2 , 
with a spectrum peaking at -240 keV (9). 

This event was extremely bright also in the 
optical and it was immediately detected with var- 
ious robotic telescopes. In particular, the Raptor 
robotic telescope detected a bright optical counter- 
part already at t = 0.5 s (10). Optical spectros- 
copy of the afterglow determined the redshift to 
be z = 0.34 (11); a UVOT UV grism spectrum 
(7) was also acquired. At this distance, the rest 
frame 1-keV to 10-MeV isotropic energy is E iso = 
8.1 x 10 53 erg, and the peak luminosity is T iso = 
2.7 x 10 53 erg s -1 . According to the luminosity 
function of Salvaterra et al. (12), we expect one 
event like GRB 130427A every >60 years. In the 
nearby universe (z < 0.4, corresponding to an 
age of -10 billion years), only a handful of long 
GRBs have been detected. These GRBs are usual- 
ly characterized by a low overall isotropic energy 
(E iso < 10 52 erg) and are associated with super- 
novae (SNe) types lb and Ic, characterized by 
broad spectral lines indicating high expansion ve- 
locities, called hypemovae (13). GRB 130427A is 
instead a powerful GRB, such as the ones typically 
detected at much higher redshifts (z > 1, with a 
mean z - 2 corresponding to an age of -3 billion 
years). The detection of a nearby and extremely 
powerful GRB gives us the opportunity to test, 
on the one hand, whether this GRB has the same 
properties of the cosmological GRBs and, on the 


other hand, whether also such bright GRBs are 
associated with SNe. Up to now, SNe have been 
associated only to under (or mildly) energetic GRBs 
in the local universe. Because a supernova asso- 
ciated with this burst, SN 2013cq, has been de- 
tected (14), we are now sure that SNe are also 
associated with very powerful GRBs, not only to 
low power bursts [fig. S6 and figure 1 of (14)]. 
Naive energetic arguments might suggest that in 
powerful GRBs, there is not enough power left 
for a strong SN; GRB 130427A definitively proves 
that this is not the case. 

The overall behavior of the x-ray afterglow 
light curve has been characterized with the main 
contribution of the XRT onboard Swift and two 
additional relevant detections from the Monitor 
of All-sky X-ray Image (MAXI) experiment (15) 
in the gap between the first and the second Swift- 
XRT observations (Fig. 2). The early light curve, 
starting from t = 260 s, is characterized by an 
initially steep decay with a slope a 0?x = 3.32 ± 
0.17, which is consistent with high-latitude emis- 
sion (16, IT); a break at t 1>x = 424 ± 8 s; and 
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Fig. 1. Details of the Swift BAT light curve in the 15- to 350-keV band. (Top) The BAT light curve 
with a binning time of 64 ms. (Inset) The BAT light curve up to 300 s, plotted on a log intensity scale, 
showing a fast rise/exponential decay feature starting at t ~ 120 s. (Bottom) The photon index values 
of a power-law model fit to the BAT spectrum in the 15 to 150 keV energy range. 



Fig. 2. Light curves for GRB 130427A in different wavebands. Swift UV and visible filters (i w2 , m2, 
wl, u, b, v); B, V, f , and /' filters correspond to Faulkes Telescope North; f , and f to the Liverpool Telescope; 
g', F, and /' to the MITSuME telescopes (7). The scaling factors for the flux density in different filters is shown 
in the inset; the scaling factor for the x-ray light curve (flux integrated in the 0.3 to 10 keV range) (7) is also 
shown. The x-ray light curve also includes two MAXI data points at t = 3257 s and t = 8821 s (empty 
squares). The fit performed over all the light curves required 24 free parameters (curve normalizations, host 
galaxy optical flux in each band, three temporal slopes, and two breaks). Because of short-term low-level 
variability superposed to the long-term behavior and possibly residual inhomogeneity of optical data taken 
from different telescopes, we added in quadrature a 9% systematic error in the optical and 5% in the x-rays. 
Final fit yielded % 2 = 543.02/565 degrees of freedom. A contribution from the host galaxy has been taken into 
account in the optical bands by fitting a constant flux of ~ 0.01 millijansky for the reddest bands, 
corresponding to r HG = 21.26 mag as tabulated in the Sloan Digital Sky Survey catalog. As a test of 
consistency, we performed the fit using a broken power law for the B, g', F, and /' filters individually. 
Uncertainties are larger, but we do find that the values of f brea k, as well as the decay indices before and 
after the break time, are still consistent, within the errors, with the values obtained by the overall fit. 


followed by a flatter decay with index ai jX = 
1.28 ± 0.01. A further break at / 2?x = 48 + 22 ks 
is statistically needed (3.8o) to account for a fur- 
ther steepening to a 2?x = 1.35 ± 0.02 (all errors 
are derived for Ay 2 = 2.7). 

The light curves in the optical and UV de- 
rived from the UVOT, as well as from ground- 
based telescopes (Liverpool telescope, Faulkes 
Telescope North, and MITSuME Telescopes) 
are shown in Fig. 2. All optical light curves are 
well fitted by a broken power law with oq = 
0.96 ± 0.01, /break = 37.4j^ 0 7 ks, and a 2 = 
1.36+° 0 ° 2 1 . Fitting the x-ray light curve together 
with the optical ones, we find the same param- 
eters from 26.6 ks onward, but to fit the early part 
of the X-ray light curve, we need another power- 
law segment with a slope of 1 .291 "q 0 02 and a break 
at 26.6+ b 6 5 ks (Fig. 2). Therefore, from 26.6 ks 
onward a common description of all the op- 
tical, UV, and x-ray behavior is possible, where- 
as at earlier times, an extra x-ray component is 
required. 

We interpret the early x-ray light curve (up to 
26.6 ks) as the superposition of a standard after- 
glow (forward shock emission) and either the 
prolonged activity of the central engine or/and 
the contribution from the reverse shock emis- 
sion (18-20). After 26.6 ks, the optical and x-ray 
light curves share the same behavior and decay 
slopes (Fig. 2), including a break at / break ~ 37 ks. 
This achromatic break is suggestive of a jet 
break, although the post-break decay (a 2 = 1.36) 
is shallower than predicted in the simplest the- 
ory [an increase in decay slope > 1 would be ex- 
pected (21)]. This could be due to additional 
components contributing to the flux, to a time 
dependence of the microphysical parameters 
governing the fraction of shock energy going 
to electrons (e e ) and magnetic field (e B ), or to 
the fact that we observed a canonical jet not 
exactly on axis, but still within the jet opening 
angle (22, 23). 

Because the optical and the x-ray emission 
belong to the same spectral power-law segment, 
it is possible to constrain the characteristic fre- 
quencies of the afterglow spectra, in turn con- 
straining the microphysical parameters of the 
relativistic shock. Additional information comes 
from the high-energy y-ray emission (2). The 
y-ray flux above 100 MeV peaks at ~20 s. If this 
emission is due to afterglow radiation, the peak 
time implies a bulk Lorentz factor r 0 ~ 500 (2, 7). 
Furthermore, the presence of the giga-electon volt 
peak suggests a homogeneous circumburst den- 
sity profile (24). Guided by these constraints in 
our choice of parameters, we used the BOXFIT 
code developed by van Eerten et al. (25) to mod- 
el the afterglow. Rather than trying to perform a 
formal fit to the data, we checked whether this 
burst, with an unprecedented data coverage and 
richness, can be interpreted in the framework of 
the standard model for the afterglow emission, 
or else whether it forces us to abandon the stan- 
dard framework. We gave more weight to the op- 
tical and higher energy fluxes because they carry 
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Fig. 3. Radio, optical, x-ray, and y-ray light curves of GRB 130427A. Radio data are from (31): 
All measurements are taken at 6.8 GHz except the later one, at 7.3 GHz. X-ray data are the same reported 
in Fig. 2 (7). LAT y-ray data are from (2). Corresponding model predictions adopt a description in terms of 
the van Eerten et ai model (25). To properly fit the radio data, only a fraction of the electrons must be 
accelerated after ~70 ks (table S10) [discussion is available in (7)]. 
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Fig. 4. Spectral Energy Distributions (SEDs) of GRB 130427A taken at different times, from 
the optical to the GeV bands. LAT data are from (2). For SED 1 and SED 2, the model is a Band 
function; for SED 3, the model is a Band + power law. Short dashed lines indicate phenomenological fits 
to the BAT + LAT data. Long dashed lines indicate results of the van Eerten et ai model (25), with 
parameters discussed in (7). The different SEDs refer to the following time intervals: 1, 0 to 6.5 s; 2, 7.5 to 
8.5 s; 3, 8.5 to 196 s; 4, 352 to 403 s; 5, 406 to 722 s; 6, 722-1830 s; 7, ~3 ks; 8, ~23 ks; 9, ~59 ks; 
10, -220 ks (host galaxy contribution has been subtracted in this case). From the optical SED analysis, the 
intrinsic extinction is negligible. 


most of the afterglow luminosity, which is orders 
of magnitude greater than the radio flux. 

Neither reverse shock nor inverse Compton 
(IC) emission are included in the model, but this 
does not affect our conclusions, which primarily 
concern the late time synchrotron afterglow (7). 
The synchrotron flux predicted by the model re- 
produces the optical emission and the x-ray light 
curve after -10 ks reasonably well, whereas the 
early x-ray flux is likely due to an additional com- 
ponent (Fig. 3). Our model underestimates the 
GeV emission, but this, given the large ratio, 
can be due to synchrotron Self-Compton emission, 
as envisaged by (2, 26, 27). The model can also 
roughly reproduce the radio emission (7). 

Consistent with the light curve analysis, the 
synchrotron flux predicted by the model reproduces 
reasonably well the optical and x-ray parts of the 
spectral energy distribution (SED) of the after- 
glow (Fig. 4), but it underestimates the GeV 
emission. Although the model does not entirely 
reproduce the complexity shown by the data, it 
does capture the main features of the emission 
properties in the pure afterglow phase. 

Overall, the properties of GRB 13 0427 A are 
similar to those of the powerful GRBs typically 
seen at z - 1 - 2 [a comparison is available in 
(7)]. This is the most powerful GRB at z < 0.9. 
It obeys the spectral energy correlations, such as 
the Fpe a k _ E iso (28) correlation and the E peak - 
Lpeak (29) correlation. Interpreting as a jet break 
the break observed at -37 ks makes GRB 130427A 
consistent with the collimation corrected energy- 
peak energy correlation (7, 30). GRB 130427A 
is also associated with a supernova, extending 
the GRB-SN connection also to such powerful 
and high-z bursts. GRB 130427A stands as an 
exceptional example indicating that a common 
engine is powering these huge explosions at all 
powers, and from the nearby to the very far, 
early universe. 
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Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial 
pulse up to 2.5 seconds is possibly the brightest well-isolated pulse observed to date. A fine time 
resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, 
consistent with models of internal synchrotron shock pulses. However, a strongly correlated 
power-law behavior is observed between the luminosity and the spectral peak energy that is 
inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the 
existing models to account for all of the observed spectral and temporal behaviors simultaneously. 

I n the context of Gamma-ray burst (GRB) Ray Space Telescope on 27 April 2013 at T 0 = 
130427A, which triggered the Gamma-Ray 07:47:06.42 UTC (2-4) is an extreme case. The 
Burst Monitor (GBM) (7) on the Fermi Gamma- peak flux on the 64-ms time scale is 1300 ± 
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1 00 photons s -1 cm -2 in the 1 0 to 1 000 keV range 
and the fluence, integrated over the same energy 
range and a total duration of -350 s, is (2.4 ± 
0.1) x 10 3 erg cm -2 . The longest continuously 
running GRB detector, Konus on the Wind space- 
craft, has been observing the entire sky for nearly 
18 years, and only one burst had a larger peak 
flux, by -30% (GRB 1 10918A) (5). GRB 130427A 
is the most fluent burst in the era starting with the 
1991 launch of the Burst and Transient Source 
Experiment (BATSE) on the Compton Gamma- 
Ray Observatory. Finally, the energy of the spec- 
tral peak in the first time bin (T 0 - 0.1 to 0.0 s), 
5400 ± 1500 keV, is the second highest ever 
recorded (6). 

The initial pulse (Fig. 1), lasting up to 2.5 s 
after the trigger, stands on its own as being so 
bright (170 ± 10 photons s -1 cm -2 peak flux for 
1 0 to 1 000 keV in the 64-ms time bin at To + 0.51 s) 
as to be ranked among the 10 brightest GBM or 
BATSE bursts (7-9). The brightness allows us to 
track the spectral evolution of the rising portion 
of a well-separated pulse with unprecedented de- 
tail (10). Evident in the GBM low-energy light 
curve [Fig. 1; as well as the 15 to 350 keV light 
curve presented in (11)] are fluctuations starting 
at around 1 s that are not present at higher ener- 
gies. If these represent additional low-energy 
pulses, their presence clearly does not dominate 
the analyses presented below. 

Past studies of time-resolved spectra of sim- 
ple pulses in GRBs indicate that there are broadly 
two classes of spectral evolution. These are called 
“hard-to-soft” and “tracking” pulses (12, 13), de- 
pending on whether the energy of the peak in the 
vF v spectrum (genetically called E veak herein) 
monotonically decays independently of the flux 
evolution or else generally follows the rise and 
fall of the flux. Typically, there are at most one or 
two spectra available for fitting during the rising 
portion of the flux history. What makes this event 
unique is that there are roughly six time bins with 
excellent counts statistics before the peak in the 
10 to 1000 keV flux. 

As seen in Fig. 1, there is a clear trend in the 
individual detector’s light curves: the >20 MeV 
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